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Abstract: We study the AC optical and Hall conductivities of Dp/Dq-br&ne intersec- 
tions in the probe approximation and use sum-rules to study various associated transport 
coefficients. We determine that the presence of massive fundamental matter, as compared 
to massless fundamental matter described holographically by a theory with no dimensional 
defects, reduces the plasma frequency. We further show that this is not the case when the 
brane intersections include defects. We discuss in detail how to implement correctly the 
regularization of retarded Green's functions so that the dispersion relations are satisfied 
and the low energy behaviour of the system is physically realistic. 



X 



* jam as@fpaxpl.usc.es 
t shock@fpaxpl.usc.es 
tarriobarreiro@uu.nl 



Contents 

[l]. Summary of results |] 

[2]. Sum rules from causal considerations ||| 

|2.1| Dispersion relations [| 

2^1 Sum rules I 



|3|. Conductivity sum rule for Dp/ D(p + 4) intersections 10 

II The D3/D7 system 



11 



]2| Dp/D{p + 4) systems [U 

Conductivity sum rule in intersections with a defect 16 

[5]. The frequency dependent hall conductivity and associated sum rules 17 

§. Concl usions and perspectives pCi 



1. Summary of results 

The interpretation of the physics encountered in the Quark-Gluon Plasma (QGP) from the 
point of view of the gauge/gravity correspondence finds an overwhelming support in the 
calculation of thermodynamic quantities and hydrodynamic transport coefficients, which 
are essential to understand the late-time behavior of the QGP in computational simulations. 
Thermodynamic quantities in the field theory are calculated by evaluating the classical 
action of the dual gravitational model on-shell. In contrast, hydrodynamic quantities 
are obtained by calculating the correlation functions, which in lattice simulations rely on 
postulating the form of the spectral density. As pointed out in Q , sum rules may be used 
as constraints on the spectral densities obtained from the gauge/gravity correspondence, 
relating the integral of -roughly- the spectral density with thermodynamic quantities and 
transport coefficients. 

In this paper we study sum rules for the current-current correlator of fundamental 
matter at finite temperature. A similar analysis has been done for M = 4SYM in 0], 
where a considerable difference between the calculations at strong and weak coupling was 
found. We will restrict to the strongly coupled regime, using the AdS/CFT correspondence 
to study different plasmas described in the string theory side by intersections of two stacks 
of Dp-branes with different dimensionalities. We will consider Nf probe Dq-bi&nes in the 
background created by N c S> Nf Dp-branes. 

The first result of the present work is the determination of the effects of having matter 
transforming in the fundamental representation on the collective modes of the plasma, 
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namely on the plasma frequency. This is related to the current-current correlator by a sum 
rule for the AC conductivity || 

Kea(oj)duj = e^LOp . (1-1) 

The AC conductivity is related to the retarded current-current correlator, G r (uj), by the 
expression ioja(uj) = G r {oj). We can then use the prescription in M to obtain the retarded 
Green's function for the holographic system and evaluate the integration of the sum rule. 
This will give us a divergent retarded correlator (and conductivity) in the large frequency 
limit. To treat this divergence we have to take a subtraction scheme in which we take the 
difference between two of these correlators. This is analogous to the Ferrell-Glover-Tinkham 
sum rule in superconductivity studies, which relates the difference between Green's func- 
tions at different temperatures (see for example ||). In this case the evaluation of the 
sum rule for the conductivity will not give us the value of the plasma frequency but the 
difference between w 2 in both setups. Schematically, we find 

■,, o, \ \-N f N e M*/\ if no defect , 

e u> 2 p (M q ,T)-e u 2 p (M q = 0,T = 0) ( x\ i q 1 , (1.2) 

[0 if defect 

with M q the dimensionful mass of the fundamental matter and A the 't Hooft coupling 1 . 
Our results show that the value of the plasma frequency gets reduced in a thermal system 
when massive matter in the fundamental of SU(N C ) is included in a brane intersection 
without defects, the reduction being proportional to the mass of the fundamental degrees 
of freedom squared. However, we are not aware of a calculation of the plasma frequency 
for a stack of Dp-branes, therefore we are calculating the effects of fundamental matter on 
an unknown quantity. It presumably scales with N 2 , whereas our correction scales with 
NfN c <C N 2 and thus is subleading in the probe approximation. 

The impact of defects 2 in the brane intersection in the system is noteworthy. Our 
results suggest that the integration of the spectral function over the frequency in the 
transverse vector channel (i.e., the integration of the real part of the conductivity) is zero 
or non-zero depending on whether the theory is a defect or not. This suggests that this is 
an example of an observable which is sensitive not only to the number of dimensions of the 
plasma in the IR, but also to the dimensionality of the UV theory (which can be reduced 
by means of a Kaluza-Klein compactification) . Though we have not been able to find a 
good reason for this difference, we have checked that this rule holds true for a large class 
of brane intersections. 

We turn our attention next to a sum rule involving the Hall conductivity. To study 
it we need to include in our setup an external magnetic field as well as a finite baryon 
density. Focusing on the D3/D7 intersection in the probe approximation, we will consider 
the Hall angle, defined as the quotient between the Hall and optical conductivities tn = 
tan## = (T xy /a xx , which gives the angle between the electric current and the electric 



1 In p 7^ 3, eo is dimensionful. Our result (^^) suggests that it is actually ~ 1/A. 

2 i.c.. theories where the fundamental matter lives only in a subspace of the original gauge theory. 
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field. Following Q we numerically determine the conductivity tensor as a function of the 
frequency and integrate tu(u) in the whole frequency range (technically we just integrate up 
to a cutoff, however the remaining contribution can be made arbitrarily small by improving 
the numerics). 

We have to deal with a holographic system which couples the modes of interest, due 
to the presence of the magnetic field, which breaks parity invariance. Going to a circular 
polarization in the transverse plane to the magnetic field we can study the system in a 
maximally coupled basis 3 which turns out to be simpler, and obtain the Hall and optical 
conductivities once the solution in the circular polarization basis is known. 

This allows us to calculate the Hall frequency, ojh, by means of another sum rule 4 

1 r°° 

UH = -P Ret H (uj)duj. (1.3) 

7T J-oo 

The Hall frequency can be understood as the frequency at which the phenomenology due 
to the Hall effect dominates over the diffusive processes in a non-ideal fluid. We present 
the numerical determination of ojh in figure |5[ At low masses of the charge carriers we 
observe that the Hall frequency is increasing, meaning that the diffusive processes become 
more important than the Hall effect for the medium frequencies of the signal. At a given 
mass value (depending on both the magnetic field and the density of charge carriers) the 
Hall frequency attains a maximum value, and for larger masses ojh decays with a power of 
the mass. 

Outline of the paper 

In section [2] we introduce some general considerations about the retarded Green's function, 
which will lead to a set of general sum rules. Concretely, we will focus on sum rules for the 
conductivity. 

In section || we investigate the consequences of the conductivity sum rule in the case of 
a Dp/D(p + 4) intersection. We reason that the satisfaction of the sum rules in these cases 
implies a reduction of the plasma frequency. We continue the discussion to more general 
brane intersections in section ||. 

In section || we turn our attention to a sum rule satisfied by the quotient between the 
Hall and optic conductivities (the Hall angle) , and we use it to estimate the Hall frequency 
as defined in the text. 

Lastly, in section |6| we write some conclusions and future directions. 



2. Sum rules from causal considerations 

The retarded correlator 

G R (t,x;t',x') = -i@(t - t') ([0(t,x),0(t' ,x')]) = -i6(t - t')p(t,x;t' ,x') , (2.1) 



3 Maximally coupled at the level of the action which leads to a completely decoupled set of equations of 
motion. 

4 Here tu = G xy /G X x is dimensionless. Thus, this equation for the Hall frequency, which is derived in 
(P . 1 9|) , yields a linear result for uj , in contrast to the case for the plasma frequency, lj p , in (1.1). 
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plays an important role in linear response theory, since it determines the expectation value 
of the operator O in the presence of an infinitesimal source <5 



(0(t, x)) a = i J di'dx'G^t, x; t', x')<5(i', x) . (2.2) 



In this paper we are interested in the case where the source takes the form of a perturbative 
electric field A u (t'). The expectation value for the current then takes the form 

(j"(*,x)) A = J dt'ax'e(t-t')([j»(t,x),r(t',x')}) A u (t',x'), (2.3) 

thus the Green's function of interest will be the retarded current current correlator. 
2.1 Dispersion relations 

The presence of the Heaviside theta function in the definition of the retarded correlator 
ensures causality in the effect of the source, avoiding signals propagating into the backward 
lightcone. For this reason the correlator naturally entering in linear response theory is the 
retarded one, since it guarantees that no signal precedes the perturbation. It also implies 
that the Fourier transformed correlator G(cj,q) satisfies the Kramers Kronig (dispersion) 
relations (see for example ]?]]) 

ReG A (u,,q)-G£ = Im ^' q W , (2-4) 

7-oo U-\L 



M*^if!^l d , (2 . 5 ) 

where V is the Cauchy principal value around the pole at [i = oj. 

The Green's function usually presents divergences in the large frequency limit, and 
therefore the results (2.4) and ( |2.5| ) may be ill-defined. This can be fixed by subtracting 



the divergent behavior in the \uj\ — > oo limit. In general one should choose a scheme 
to perform this regularization. One usual choice in holographic contexts is to take off 
the zero temperature Green's function, given by powers and logarithms of the frequency 
(Ga = Gt — Gt=q)- It may be that the exact expression for Gt=o is n °t known, in which 
case the comparison should be carried out with a specific scheme, and the interpretation 
of the results has to take this scheme into consideration. It turns out that for physical 
quantities we can show that there is no scheme dependence. 



Therefore, in order to apply ( |2.4| ) and (2.5), we have to make sure that we have 
subtracted all the divergent contribution to the correlator at large frequencies [|l], |8|. We 
have marked this in the previous equation by the a subscript. 

It may be that although the subtracted function Ga has had all UV divergences re- 
moved, a constant piece may still survive (ie. Ga(w — > oo) 7^ 0) jl] and thus the integral 
in the Kramers Kronig relation diverges. This constant (real) piece, G^, which remains in 
the UV, corresponds to a contact term given by a 5(t) factor in the two-point correlator 
function, and is not strictly acausal. This factor may be obtained by OPE techniques, re- 
lating it to the expectation values of local operators The Kramers Kronig relation then 
holds for the fully subtracted Green's function which has had the contact term removed. 
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Causality of the Green's function, as dictated by the Kramers Kronig relation, de- 
termines the symmetry properties of the spectral function. In general, the spectral func- 
tion of bosonic hermitian operators p(ui, q) is an odd function of the frequency defined 
as pa = — 2ImGA, where pa stands for the part of the spectral function that survives 
after the subtraction (NB. The imaginary part of the Green's function will not receive any 
contact term contributions, which effect solely the real part). 

2.2 Sum rules 

It is possible to extract a number of relations between integrals in the whole frequency 



regime and some physically relevant quantities. These relations follow from (2.4) and 
and involve on the one hand an integration over the entire positive and negative frequency 
domain, capturing the effects occurring at infinitely small (positive and negative) times 
and on the other, they imply a certain behavior over the whole frequency range, including 
for null frequency, thus giving information about the response of the system at large times. 

Recall that we have already subtracted the large frequency behavior in the integrand 
to make certain that these integrals converge. A stricter method would be to introduce a 
cutoff frequency, uj\, that effectively neglects instantaneous effects in the response of the 
system. In fact, in this paper we will perform the integrations numerically, such that this 
cutoff frequency appears naturally in the calculation. Nevertheless, the dependence of the 
results on the position of the cutoff (once we have subtracted the non-convergent part) gets 
suppressed as ui\ — > oo, and the difference between cutting the integrals off or performing 
the full integration becomes negligible. 

The information we gain at null frequency (i.e., infinitely large times) has to be in- 
terpreted with some care. Here, we will restrict our studies to the probe approximation, 
in which we insert a probe in a never-changing background. Physically, neglecting the 
backreaction of the probe on the background is only valid at times lower than a specific 
timescale, after which the effects of the probe on the system must be taken into account 
(i.e., the time at which the transfer of momentum from the NfN c degrees of freedom to 
the iV c 2 de grees of freedom becomes important). Therefore, we should interpret our results 
as valid for time periods well beyond when the perturbation is turned on, but not as large 
as for the probe approximation to fail. 

Thermodynamic sum rule 

Perhaps the simplest sum rule we may obtain from the dispersion relations at zero momen- 
tum 5 is the value of the retarded correlator at zero frequency. Inserting p~A = — 2ImG , A in 
dU) we obtain 

ReG A (0,0)-G% = -^rr PA ^ ,0) dw, (2.6) 

which is a thermodynamic sum rule for the spectral function. This is equation (11) in Q, 
with a relative factor of 2ir in the definition of the spectral function, and a global sign 
coming from the definition of the Green's function (see also |TJ). To obtain this relation 



'For the rest of this paper we restrict to the q = case. 
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we used only the general definition ( |2,l[) , so it does not depend on any details of the model 
under consideration. 

In Q certain sum rules for the spectral density in hot gauge theories were derived. 
Specifically, for the M = 4 SYM at infinite 't Hooft coupling we have the following sum 
rules involving the xy, xy graviton correlator (following our conventions, which differ from 
the ones in |l| by a factor of 2 in the definition of the spectral function) 

1 />oo J 

Re G™ (0,0) - G™ = — / —p x ^ xy (u) , (2.7) 

^ Jo 10 

where Re G x ^ ,xy (0, 0) = P is the pressure (from hydrodynamics, see H) and G^ = ^-P is 
a constant value that survives at infinitely large frequencies |IJ. Subtracting and using the 



conformal equation of state, the l.h.s. in (2.7) is — f e, where e is the energy density. 



In |2| a sum rule was obtained for the R-current correlator at vanishing momentum in 
TV = 4 SYM at infinite 't Hooft coupling 

i poo ~R— charge / \ 

Re Gl- char9e (0) - G™ = — / ^ [ ^-doj = , (2.8) 



7T 







where = and ReG^ char 9 e {^ = q is a result dictated by hydrodynamics. Other 



examples of such sum rules can be found in [11, 12]. 



Recently, the authors of |T(| generalized these sum rules by studying the constraints 
of gravity on the retarded correlator. Our interest in this paper is to use the sum rules 
to study the behaviour of transport coefficients, since the systems we study are known to 
have a meromorphic retarded Green's function with causal structure. 

Conductivity sum rule 

From ( |2.3| ), and using Ohm's law, we can express the AC conductivity in terms of the 
subtracted retarded correlator as 

Ga(w) = iuja(u;) . (2.9) 
Therefore, the pole structure of the conductivity is basically the same as in the retarded 



correlator, and the dispersion relations (2.4) and (2.5) have to hold for a(u). There is the 
possibility of an additional pole (accompanied by a delta function) at zero frequency in the 
conductivity. However, this pole does not affect the satisfaction of the dispersion relations 
discussed above. Consistency implies that the cutoff-dependent integrated conductivity 

W(lo a ) = - / a(n)dn, (2.10) 

n J-WA 

goes to — G^ in the limit wa — > oo. To see this, notice that from ( ^?9| ) Re G R {ui) = 
— uolm a(uj) and lmG R (uj) = u:Rea(uj). Plugging this into (p.4|) we obtain 



-clmaM -Gt = ~-Vr ^^W , (2.H) 

J-oo W-jU 
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and from ( |2.5D for a 



-wlnifrw = P/ ^ du. (2.12) 

7T U-fi 



As the imaginary part of the correlator goes to zero at large frequencies, i.e., G^ is 
real, then a°° = follows. Assuming now that in the principal value all the limits are 
taken at the same rate, we can join the former integrands to obtain 6 

-G& = -- vT Rea W / _ w)d ^ = I f°° Re(T ( M )d/i . (2.13) 

Now, as Re G R is even in frequencies, the corresponding quantity Imcr will be odd, meaning 
that the integrated conductivity ( p. 10 ) is real. Therefore 



lim W(u A ) = -G%. (2.14) 

The integral over the whole range of frequencies is employed in the analysis of plasmas as 
a way to obtain the plasma frequency, giving the energy of the collective excitations in the 
plasma, by using the expression H| 

lim W(uk) =e ul. (2.15) 

Because in our case we are defining the conductivity from a subtracted correlator, what we 
are calculating is a specific contribution to the plasma frequency (in fact it is a difference 
in contributions as we shall discuss shortly). Notice that we are not demonstrating here 
that the integral of the conductivity in the frequency domain determines the value of the 
plasma frequency. We are using this as a definition of the plasma frequency to calculate it 
in a holographic setup, instead. 

It seems important at this point to summarize our results. The causal structure of 
the conductivity allowed us to write down equation fl2.4|) . For co = this relation gives 
us the difference between the Green's function at infinite frequency and at zero frequency. 
On top of this we use the fact that the Green's function at zero must vanish so that the 
real part of the conductivity does not develop a delta function at the origin. A delta 
function in Re\a{uj)] = it5(uj) is not a sign of superconductivity, but is a universal property 
of systems with translational invariance (see for example jOJ). In systems where the 
quenched flavors are modelled with probes, the huge mismatch in the number of degrees 
of freedom N c > Nf mimics the presence of fixed scattering centers in the bulk, leading 
to a finite DC conductivity. In reality this is a transient effect for times of order Nt/N c 



(see [14]). Of course, in a setup where backreaction of flavor is taken into account, full 



translational invariance is recovered, and a delta function should be present. It would be 



worth studying this effect, for example, in the recently constructed model in [15|. For us 



6 Note that plugging the relation between Re a and Pa(w) in (2.13) we do not recover (2.6) with 



ReGA(O) = because the Cauchy principal value to evaluate the pole at ui = has disappeared. In 



particular, the r.h.s. of equation (2.6) is insensitive to delta functions 
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here, the absence of a delta function is a physical constraint on top of causality. This 
definition then gives us a way to calculate not only the difference between the zero and 
infinite frequency behaviours, but the absolute value of the Green's function at infinity. 
This itself has a physical interpretation as the plasma frequency. In a metal the plasma 
frequency describes the frequency of vibration of electrons about the ions, however in our 
system there is clearly nothing analogous to the ions. A possible interpretation of the 
plasma frequency may then be given by the rate of pair production by an infinitesimal 
electric field which causes pairs of charged particles to be pulled out of the vacuum before 
annihilating again. The fact that the plasma frequency is dependent only on the quark 
mass and not the baryon density suggests that it's a vacuum effect and not related to the 
excess of charge carriers. 



Expression ( 2.10 ) relies on the relation (2.9) and the causal structure of the conduc- 



tivity (and thus the retarded correlator). See however [16| for a discussion deriving similar 
relations from the gravitational point of view. Equation (2.10) can also be used to ob- 
tain numerically the weight of the factors multiplying delta functions in the conductivities. 
These delta functions can occur in metals not presenting superconducting behavior, but 
that are not able to dissipate momentum. In this case, impurities will destroy the delta 
function giving rise to a Drude peak. In superconducting systems this is not the case, the 
delta functions will still be present in the presence of impurities, but their weight may get 
reduced by the presence of the extra scattering centers. 

Hall angle sum rule 

When turning on a magnetic field the Green's function of the transverse vector channel 
has to be treated as a matrix, since the parity-breaking associated to the finite magnetic 
field induces a coupling between the fluctuations along the transverse directions to it. In 
this case, the Kramers Kronig relations are written as 



Wf.A nA ( ,.,_ i r Gl{n)-Gl 



G%(u)-G% = -I ^^dfi, G£(u) = -I " A d/x, (2.16) 

where the superindices ' denote hermitian and antihermitian parts, respectively. In this 
case G^ is a constant hermitian matrix. 

As we will see later, the presence of a magnetic field in the homogeneous system 
induces a form for the Green's function matrix (notice that now we are not subtraction the 
asymptotic behaviour) 

G = (% x G jA • (2.17) 

\ ^xy ^xx I 

We will focus our interest on a sum rule satisfied by the quotient tff(w) = G xy /G xx . Let 
us sketch the procedure following {?], § . 

We would like to show that tn is a causal signal. If this is the case then tn cannot 
have poles in the upper half plane. Now, the x — y component of the Green's function 
matrix has to be a causal function too, therefore potential acausal poles in tjj must come 
from zeros of the diagonal component in the upper half plane. However, G xx {uj) is a 
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meromorphic function constrained to satisfy G xx (ui)* = G xx (—uj*). In this case, complex 
variable analysis leads to the conclusion that G xx (u) takes real values only on the imaginary 
axis, where it behaves monotonically from G xx (0) to G xx (zoo) (see 0, chapter XII). As we 
will show, G xx (0) = in the cases of study and at infinite frequency it diverges, therefore 
the only pole coming from the zeroes of G xx is at the origin and can be avoided by choosing 
an appropriate integration contour in the complex plane. 

Now, following ||, consider the Kramers Kronig relation 7 

, , v, 1 f°° uiRe\t H (a)} 
uIm[t H (u)] = -V 1 "^ d/j. (2.18) 

7T J-oo U- H 

Defining ujh = lim^^oo ui Im[^(w)], in the large uj limit we obtain 

1 r°° 

uh = -V Re[t„ (2.19) 

7T J-oo 

This constant ujh is interpreted as the Hall frequency in || and gives the rate of the preces- 
sion of charged particles in the magnetic field (note that with the AC current the charged 
particles will only move around a small part of this circle during any given oscillation of the 
electric field). Thus the Hall angle together with the Hall frequency can tell us about the 
helical motion of a particle in the setup of interest as discussed in more detail below. As 
in the case of the plasma frequency, we use the sum rule to define the Hall frequency and 
calculate it in a holographic setup; we are not demonstrating that both quantities are the 
same. This would have to be checked with an independent calculation. It is also possible 
to apply the hydrodynamic sum rule to this quantity, obtaining 



ReM 0)] = a -BL 



w=0 



-i r-oo j 

-V — Im[iffH], (2.20) 

vr J-oo u 



where the DC ohmic and Hall conductivities can be calculated analytically in most cases. 

We can understand the definition for the Hall angle by relating it to an injected current 
in the medium along, say, the direction x 

_ G xy (u) a xy {u) E x j y (u) 

tH(U) - ~ - = / \ T-i — i—r , (2.21) 

G xx (u) a xx (uj) E x j x (uj) 

where j y is the induced Hall current along the direction y. Fourier transforming back into 
time domain we find 

/oo 
t H (t - t')j x (t')dt' , (2.22) 
-oo 

with tnit) the Fourier transformed Hall angle. Therefore, the presence of a magnetic field 
forces the electric current to precess in a dynamic process described by the Hall angle. The 
hall frequency ujh determines the value of \m\ t _± t i ti{(t — t'), hence describing this precession 
at very short times, and for a perfect fluid, where diffusive centers are absent, should be 
equivalent to the cyclotron frequency || . In the case we study in section || this is not the 



7 Since tu = G xy /G xx and in the large frequency limit G xy — > and G xx diverges, then — 0. 
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case though, since we consider a probe approximation and the fixed background acts as an 
infinitely massive diffusive center. 

As mentioned before, tn is analytic in the upper half-plane. The integral will converge 
in the case we study later since the real part of tn decays faster than X/oj at large frequen- 
cies 8 . However, this will not be true for the imaginary part, given the finiteness of the Hall 
angle from figure || and the definition of the Hall frequency ujh = hnioj^oo oj Im[i#(u;)]. 

3. Conductivity sum rule for Dp/ D(p + 4) intersections 

The metric and dilaton corresponding to the decoupling limit of a stack of N c S> 1 black 
Dp-branes is given in the string frame by [17] 



ds 2 = H -l/2 ( _ /dt 2 + df 2) + H l/2 /dT^ + r 2 dQ 2_\ e <P = ff^ , {3A) 

where / = 1 — (rh/r) 7 ~ p and H = (L/r) 7_p . The additional RR form does not play a 
role in this paper. The matter of the theory transforms in the adjoint representation. By 
considering an intersection of the stack of N c .Dp-branes with a stack of Nf <C N c probe 



-Dg-branes matter transforming in the fundamental representation can be added [18|. It 
was shown that, in order to preserve supersymmetry in the extremal case, —> 0, the 
dimensionality of the fundamental branes, q, has to be given by q = {p, p + 2, p + 4}. In 
the latter case the probe branes wrap the worldvolume directions of the Dp-branes plus 
an (S 3 contained in the S 8 ~ p appearing in ( |3.1| ) whilst in the other cases the intersection 
must contain a defect in the worldvolume directions of the Dp-branes. Working in a 
dimensionless radial coordinate u = r^/r the intersection in the latter case is summarized 
in the following table: 



t W u S 3 X, 



5—p 



Dp x 
D(p + 4) x x 

with the dynamics governed by the DBI action 

Sdbi = -N f T D(p+A) [ e-V-det [5 + 271-^], (3-2) 

JMp+4 

where T D ( p+4 ^ = 1/ ( y {2TT£ s ) p+i g s l s s ) is the tension of the brane, g the pull-back of the 10- 
dimensional metric to the worldvolume of the probe brane .M p +4, t s the string length and 
F the field strength of the gauge field living on the probe brane. 

The embedding profile of the probe branes can in general be specified by a single scalar 
function ip{u) which is related to one of the angles in the compact X^- p manifold (there 
are situations where the embedding must be written in terms of multiple scalar fields but 
we here specialize to the more generic cases). This can be specified by decomposing the 
metric of the (8 — p)-sphere as 



dQ s-P = + + (1 " • (3.3) 



s For this, it is crucial that we do not subtract the large frequency behaviour of the optical conductivity. 
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In this case, the asymptotic behavior of the embedding function near the boundary, u — > 0, 
is given by ip ~ m q u + c q u 3 . The dimensionless constant m q is directly related to the 
quark mass (see ( 3.1 2| ) ) and c q to the quark bilinear condensate |ll 



The system presents two different phases, separated by a first order phase transition 
at a critical mass m* [jnj. For masses below the critical value a horizon is induced in the 
worldvolume of the flavor branes, giving rise to a conducting phase in which the mesons 
become unstable. For masses larger than the critical value the transverse S 3 wrapped by 
the branes shrinks to zero size above the horizon of the background geometry, and the 
system is in an insulator phase, with a stable mesonic bound-state spectrum. 

It is possible to introduce a finite baryon density into the description of the field theory 
by turning on the temporal component of the U(l) abelian center living on the worldvolume 
of the probe branes. On the horizon this component has to satisfy Aq(1) = 0, in order 
to avoid a singular one- form. From the asymptotic behavior on the UV, Aq ~ fi — au 2 , 
we identify \x as the chemical potential and a as related to the baryon density. Whenever 
a / 0, the embedding has to be of the black hole kind, but a phase transition still remains 



if the value of the baryon density is low enough [ 20 1 . 

In the present work we are interested in the fluctuations of the U(l) vector field in 
the conducting phase. The expression (|2.6| ) is evaluated at null spatial momentum for the 
fluctuations, so an SO(p) symmetry is preserved and the p spatial components of the vector 
field can be studied by considering just one of them, say 8A p (t, u) = exp(— iut)5a p (u), with 
5A the perturbation of the gauge field. We will work in the 5a u = gauge. The equation 
of motion for the fluctuation is obtained by expanding ( |3^ ) to quadratic order, and is 
integrated imposing that it is an ingoing-wave at the horizon of the black hole and fixing 
an overall normalization. In the UV it behaves as a p (u ~ 0) ~ A p + B p u 2 + • • • . An 
additional logarithmic term can appear in this expansion. We will see later how to deal 
with these. 

One can now obtain the two-point function & > J (u, 0) using the standard prescription 
(from now on we will focus on the J p -J p correlator, so we will skip this superscript) 



(3.4) 



where V(u) is proportional to the wronskian for a p , and can be read from the (d u Sa p ) 2 
term in the quadratic action. Notice that it captures several constants like 2V/Td Cjj+4 . , 
(2tt£1) 2 or the volume of the S 3 . 

3.1 The D3/D7 system 

We first specialize to the D3/D7 brane intersection, corresponding to quenched fundamen- 
tal matter added to four dimensional (thermal) M = 4 SYM with a small number of N = 2 
hypermultiplets. 

In order to evaluate ( [2.6| ) we must subtract the T = contribution from the 'bare' 
Green's function. We subtract the analytic part of the correlator at zero mass and zero 
baryon density which can be found in pi]]. Evaluating the retarded correlator at large real 
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frequencies we obtain 



G s {u) = --^-(^-\ \2 lE + \og [e 2 l-zvrsignH]. (3.5) 



N f N r T 2 ( 

,2vrry i " 1 ~° \~ V 2?rT , 

The e — >• limit should also be understood in the previous expression. The logarithmic 
divergence is avoided by the use of holographic renormalization. To make this more concise 
we can Fourier transform ( |3.5| ) back into the time-domain, obtaining a correlator 

Gs(t) = ^ Q log(e 2 )^) + nd 2 ^j , (3.6) 

and we see that the T = contribution is a retarded correlator with a contact term, 
which presents the logarithmic divergence. The 5(t) contact term could alternatively have 
been removed by adding the appropriate counterterm to the DBI action via holographic 
renormalization given by 



Lp = -L N Td7 V(S 3 )(2t:£ 2 s ) 2 [ V^h^h^FijF^logie 2 ) , (3.7) 
4 JMa u ^ e 

where is the 4-dimensional metric restricted to the u = e slice. The subtraction of 



( |3,5| ) is equivalent to the addition of this counterterm. However, the term proportional to 
d 2 @(t)/t cannot be removed with holographic renormalization and is needed to regularize 
the spectral function at large frequencies. 

Having tamed the divergences, we can investigate the finite temperature correlator as 
a function of the frequency in a representative range of masses and baryon densities 9 . We 
present a typical profile for the retarded correlator in figure [l], and we will refer to this 
figure as representative of the whole set of correlators we have studied. 

One first thing to notice in figure [l|(left) is that ReGA(O) = 0. This is a result 
that has been obtained analytically in |23|] to be G(oj) = iuhjdc + 0{uj 2 ) } just reflecting 
the consistency of the numerical calculation of the retarded Green's function with the 
hydrodynamic limit. In the former expression gdc is t ne conductivity of the system, 



which coincides with the macroscopic calculation in |22]. 

We need now to calculate the value of the retarded correlator in the infinite frequency 
limit and check that it coincides with the r.h.s. of (|2.6[). Studying this for different values 



of the quark mass and baryon density we find that the results are consistent , and that 



Gt = ReG a( 0,0, + I„f £^«^> d . _ M*L K _ (3 . 8) 
7T /_„ ui 6 H 3A H 



with M q = hvATrriq, the dimensionful mass p^]. This result, up to a factor 2, is the same 



as that obtained in appendix B of \21\ for the T = scalar-scalar correlator. In this limit 



both the scalar and vector channels are related by super symmetry. 



9 



The dimensionless baryon density d that we specify in the examples is related to the a constant in the 



2 ~ 



asymptotic expression of the Aq mode by a — d. 

10 We found differences between the numerical and analytic values of only O(10 -3 — 1CP 4 ) 
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Figure 1: Real (left) and imaginary (right) parts of the correlator function for finite mass m q = 
1.3266, finite baryon density d — 2 and zero momentum, in units of N c NfT 2 /2. The dashed line 
marks the value rn^/3 « 0.58662. 



One could argue now that our subtraction scheme is incomplete, and that we should 
have also included the part in (|3.5|), which corresponds to a contact term proportional 
to 5(t). This would have changed the value of RcGa(O) to a negative value, a fact with 
two important consequences. The first one is that the result would have been inconsequent 
with the hydrodynamic calculation, since this calculation is insensitive to the existence of 
such a contact term. The second consequence is however more crucial. A negative value 
of the correlator at zero frequency reflects, via the relation with the conductivity, in a 
1/uj pole in the imaginary part of the conductivity. The Kramers Kronig relations tell us 
that we should therefore find a delta function at zero frequency in the real part of the 
conductivity. This situation, however, cannot be realized within the probe approximation, 
since the delta function signals a saturation in the energy absorption of the system, but in 
the strict probe approximation limit the background is always able to absorb the energy 
pumped by the fields on the probe branes. 

Coming back to ( |3.8| ), in the m q — > limit we recover the R-current result 11 given 
by equation Q2.8| ). Also, this expression holds independent of the quark condensate or the 
baryon density. In particular, this means that the relation does not feel the unstable phases 
present in the vicinity of the phase transition, and will still hold if we overheat the system 
beyond the critical point, provided we do not reach the temperature at which the first 



quasinormal mode becomes tachyonic [19|, since then the dispersion relations will not hold 



due to the presence of poles in the upper half plane. 



We have also checked that the relations ( |2.4D and (2J3) hold, obtaining numerically 
figure |l](left) from |l](right) and vice versa. This can be used as a check of the numerical 
calculation, since it means that the signals in the system propagate causally, a fact that 
numeric instabilities in the integration of the equation of motion could spoil. 

An example of how numerical instabilities do not affect our results to any great extent 
is given by figure ^. In this figure we show the absolute value of the integrand in (| 



1 In this limit the R-current results are obtained by the rescaling NfN c — ¥ N% 
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2 n | p A (u») | 
N c N f T 2 vo 

0.0773 -1.2395 




1 2 3 4 5 

Figure 2: Logarithmic plot of the integrand for the sum rule (in units of N f N c T 2 / (2tt)) as a 
function of to = uj/(2ttT) for m q = 1 and d = 0.1. The amplitude of the oscillations decays as 
ijj 2 e~ kul . The numbers above the line are the contribution to the integral by each oscillation, adding 
up to approximately — ir/3, in agreement with the dispersion relation. 



for a given value of the quark mass and baryon density, together with the value of the 
integration in the positive frequency axis. Comparing the value of this integral with the 
value given by (|3.8[) we obtain an overwhelming agreement, meaning that although we are 
not numerically integrating to infinite frequencies, the vanishing error introduced does not 
affect our results or conclusions. This is because the amplitude of the oscillations around 



the subtracted T = contribution decays exponentially, as observed in [24, || 

In figure ^ it is seen that the last plotted bump already presents some small numerical 
errors. This error can be decreased with improved numerical methods. 

Integrated conductivity 

The results described above can be expressed in terms of the sum rule for the conductivity 
(2.10) by stating that 



W mq = lim W(u A ) = - NfN (i cT2 m 2 q , (3.9) 

where we have specified the value of the mass of the fundamental matter. Recalling the 
sum rule djJD we observe that the effect of introducing massive fundamental matter in the 



system is to effectively reduce the plasma frequency. The more massive the fundamental 
matter, the bigger reduction we get, though clearly this is a subleading effect in Nf/N c . 

3.2 Dp/D(p + 4) systems 

We will restrict the discussion in the current section to p = {1, 2, 3, 4}, since the back- 
ground created by Dp-branes with p > 5 is thermodynamically unstable |25| and we cannot 
define the current-current correlator at zero momentum, equation ( |3.4D , for p < 1. These 
setups are the holographic duals to p+ 1-dimensional theories, and therefore its interest is 
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mostly academic in the D4/D8 case (note that we are not introducing compact directions 
for the Dp-brane world volume directions), p = 1 and 2 may be used to model some 
condensed matter systems. 

We have not found in the literature the equivalent to the asymptotic expression ( |3.5| ) 
for p = {2, 4}, which is needed to define the subtracted Green's function Ga(w) entering 
in the expression for the sum rules. When p = 1 the integration of the vector mode can 
be performed analytically in the massless, zero baryon density case, obtaining a solution 
satisfying the ingoing-wave condition at the horizon and normalized at the boundary 



ai(it) = exp 



3a; f u udu 

-i- 



2n T /n 1 - ifi 



(3.10) 



From here it is straightforward to see that the retarded correlator is G(uj) oc —i-^pf- This 
is precisely the quantity we should subtract, since it diverges at large frequencies and 
corresponds to nothing but a 8(t) contact term. It is then clear that, in this case, the 
sum rule in the massless, zero baryon density case is zero, since the integrand vanishes 
identically. We found the result for the massless sum rule to be independent of the baryon 
density. 

We do not have an equivalent expression to (HD or ( glcj ) for the D2/D6 and D4/D8 
cases. Despite of this, we can evaluate the conductivity sum rule by using experience 
learnt from the D1/D5 and D3/D7 systems. This we can do by changing to an alter- 
native subtraction scheme, in which we do not only subtract the analytic divergent part 
of the correlator at T = 0, but take away the whole Green's function corresponding to 
the massless, zero baryon density case at finite temperature. With this, we are adding 
some frequency-dependent structure to the spectral density of Ga, originating from the 
introduction of an infinite set of quasinormal modes (in the complex plane). 

The new subtraction scheme redefines Ga(uj) = G(u;,m q ,d,T) — G(u;, 0, 0, T), as 
opposed to the definition given in section \2.l[ which assumed also T = in the sub- 
tracted piece. Therefore, the conductivity sum rule ( 2.10| ) evaluates now to W m j = 



—G^(m q , d) + G^(0, 0). However, if the aim is to compare the integrated conductivity for 
two different masses, this shift will have no effect and our results are not affected by choos- 
ing this subtraction scheme. In the D3/D7 and D1/D5 cases we found that (0,0) = 0, 
and we will consider that this holds in the following. This seems a reasonable assumption, 
since the system has not a natural scale to provide the dimensions for G^O, 0). 

Performing a scan of the mass and baryon density parameter space, we find that 
the large frequency behavior of the integrand in ( p. 10 ) again decays exponentially in the 



D4/D8 case, whereas in the D1/D5 and D2/D6 systems it decays just with a power 
of the frequency (ui~ 2 and w~ 7 / 3 respectively). However, these systems are more stable 
numerically, and we can obtain reliable results at large frequencies. The calculation of the 
spectral sum rule can be summarized by the expression 

W mq = -2irTaoK p m 2 q , (3.11) 

where do = NfN c M p is the DC conductivity at zero mass and zero baryon density, and 
the values for M p and k p are given by 
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p 




K 


1 


0.598(3) 


1/VttA 


2 


0.713(2) 


(9T/160vrA) 1/3 


3 


0.333(1) 


T/4vr 


4 


1.199(1) 


AT 3 /54vr 



Hp is the parameter we have determined numerically. Notice that the conductivity 
has energy dimension 5 a = p — 2, whereas the 't Hooft coupling, A = N c gy M , has energy 
dimension 5\ = 3 — p in each case. These expressions can be written in terms of the 



dimensionful mass [19] 



such that 



M n 



W Ma 



Hi 



9-p 

27-p nl 



-m 



Q J 



Ml 



-PpN f N c — , 
2k 3 and fa = 3k 4 /2 2 / 3 . 



(3.12) 



(3.13) 



where p 1 = Ski/2 1 / 3 , fo = 5k 2 /2 6 / 5 , f3 3 
4. Conductivity sum rule in intersections with a defect 

Up to now we have studied theories where the intersecting dimensions of the Dp/Dq system 
were precisely those of the L>p-brane worldvolume, however, it is possible to study similar 
setups, where a defect does exist in the intersection. One of the possibilities is the Sakai- 



Sugimoto model [26|, which has received lots of attention as a holographic model of QCD. 
However, it is not clear how a finite quark mass would be added, so we cannot use this 
system to find a similar result to (|3.11 ). Evaluating the sum rule in the massless case we 
obtain that the integrated conductivity vanishes. 

Another candidate system to describe properties of 3 + 1 plasmas is a D4/D6 intersec- 



tion with one compactified dimension [19]. This intersection is summarized in the following 
table 





t 


R 3 x 4 


u S 2 X 2 




X 


X X 




DQ 


X 


X 


X X 



As in the D3/D7 case, it is possible to add a finite baryon density and perturb the gauge 
field to obtain the retarded current-current correlator. We have performed a study calcu- 
lating the sum rule for different values of the mass and the baryon density and found that 
the integrated conductivity vanished always. This suggest that the presence of a defect 
dimension in the theory has a very important effect on the determination of the sum rules. 

We studied also the case of a 2 + 1-dimensional theory with one defect, namely the 
D3 1 D5 intersection with a non trivial flux along the defect to stabilize the setup |27]] . The 
intersection is 





t 


R 2 ^ 3 


u S 2 X 3 


D3 


X 


X X 




D5 


X 


X 


X X 
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Considering again the current-current correlator and scanning the parameter space of the 
theory for several values of the mass and baryon density we find that the integrated con- 
ductivity always vanished in contrast to the non-defect theory. 

Hence, we find that in flavor models based on intersections at a defect dimension there 
is no shift in the plasma frequency due to the presence of fundamental degrees of freedom. 



5. The frequency dependent hall conductivity and associated sum rules 

In this section we calculate the frequency dependent Hall conductivity and associated sum 
rules in the case of the D3/D7-hvane intersection 12 . In order to find the Hall conductivity 
in a macroscopic setup one must turn on both an electric and perpendicular magnetic field 
and find the current induced in the transverse direction to these external fields. Using 
Ohm's law one may then deduce the conductivity. However, such a setup will in general 



only provide the DC hall conductivity. This has been done previously in [29]. Here we are 
interested in the frequency dependent (AC) conductivity, for which a microscopic calcula- 
tion is needed. The setup is therefore as follows. A probe Z)7-brane is embedded in a black 
D3-brane background and a gauge field of the following form is turned on 

A = ^^-dt + e- iwt 5a x (u)dx + ( e~ iu]t 5aJu) + ^-^5 x] dy . (5.1) 
2ira' \ zira' J 

The gauge field therefore encodes a magnetic field in the z-direction, a baryon density 
given by the temporal component of the gauge field and excitations of the gauge field in the 
x and y directions. Note that the choice to include only zero momentum modes will simplify 
the calculation a great deal as it allows us to decouple scalar and vector fluctuations. 

The combination of the magnetic field and the finite baryon density causes a coupling 
between the x and y components of 5a^, and thus solving the differential equations be- 
comes much more complicated. However, it is possible to chose an appropriate basis for 
the gauge field in which we can decouple the equations of motion. This is the basis of 
circular polarizations which we label ai(u) and an(u), standing for left and right circular 
polarizations. This change of basis is implemented by 

5a L (u) + 5a R (u) 5a L (u) - 5a R (u) 
5a x (u) = , 5a y (u) = . (5.2) 

From the usual DBI action, we Legendre transform with respect to the time component 
of the gauge field at zeroth order in the fluctuations in order to calculate the embedding 
profile for the D7-brane. Such a calculation has been shown in detail numerous times 
(see for example J20[]). In all of what follows we will be at finite baryon density, as the 
Hall conductivity is only present in this case. In the presence of finite baryon density all 
embeddings are of black hole type and end with the Z)7-brane on the black hole horizon. 

In order to calculate the fluctuation solutions for the L and R components of the gauge 
field, we expand the action to quadratic order in these fields allowing a plane wave ansatz 
for the solution. Off-shell there are terms in the action which are quadratic in both the 



12 



See |^8| for a similar calculation in a setup with a defect. 



-17- 



5dL and 5a r fields as well as cross terms which mix the two. However at the level of the 
equations of motion, where the embedding solution is put on-shell, this coupling in the 
equations vanishes (i.e. the coupling terms in the fluctuation e.o.m. vanish if we plug the 
solution of the embedding profile). In fact the coupling at the level of the action becomes 
pure and the L component equation of motion involves only the R fluctuation and vice 
versa. 

Having uncoupled the equations of motion we can calculate the Green's function with a 
numerical integration, by specifying the incoming wave boundary conditions at the horizon 
and solving out to the UV where we read off the asymptotic behaviour. However, although 
the fluctuations decouple in the equations of motion, they couple exactly (there are no 
non-coupled terms) in the on-shell action. The Green's function, as a matrix in the space 
of left and right circular polarizations, is therefore exactly off-diagonal. The UV behaviour 
of the solutions is of the form 

5a L ,R( u )\u^o = ai,R + b L ,RU 2 . (5.3) 

By studying the UV behaviour of the action therefore one can determine the exact form 
of the Green's function matrix, which is given by 



Gll Glr 
Grl Grr 




(5.4) 



/ G X x 


G X y j 


-i( 


\G yx 







(5.5) 



We can study the matrix of Green's functions in the original (x, y) basis by transforming 
the matrix in the (L, R) basis, which is calculated from the numerical calculation. The 
(x, y) Green's function matrix is then given by 

Grl + Glr "i{Glr — Grl)\ 
-i(Glr — Grl) Grl + Glr J 

This defines the bare green's function, and so in order to study the Kramers Kronig relations 
we must subtract the large oj divergent behavior as discussed in the previous sections. 
Notice that, naturally, G xx = G yy and G xy = —G yx . 

The equations of motion for Scxr and Sol 13 are equal up to a sign change of B, as 
expected, and thus the above Green's function is governed by a single equation, where 
we must chose the appropriate sign of B to perform the calculation. From the Green's 
function the conductivity follows simply and thus we are left with the final expression for 
the complex hall conductivity 

&xy = "j— • (5.6) 

Aoj \a L a R J 

Note that the large oj divergence cancels in this expression and thus the large oj subtraction 
is unnecessary. This we will find not to be true for a xx ■ We can write the Hall conductivity 
alternatively as 



1 I 5a' L (u,B,d) 5a' L (u,—B,d) 



Auoj \5a L {u,B,d) 5a L (u,-B,d) 



(5.7) 

u->0 



13 Details of the equations of motion can be found in equation (A. 2) in 
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Figure 3: (Left) Real part of the regularized Green's function for fixed mass m q — 0.297, zero 
baryon density and increasing magnetic field B = (blue), B — 2 (red) and B — 5 (yellow). 
(Right) Real part of the regularized Green's function for fixed magnetic field B = 5, zero baryon 
density and increasing angle of incidence of the embedding profile in the black hole, corresponding 
to masses m q = 0.297 (yellow) and m q = 0.256 (green). Note that the values of tpo are 0.72 and 
0.9 respectively. For the larger mass, smaller ipo, the Zeeman effect is hidden by the broader peaks 
in the spectral function. 



where the primes are derivatives with respect to u. There are some immediate points to 
note about the calculation and the preliminary results. As explained in previous sections, 
the Green's functions should satisfy a Kramers Kronig relation. This is both a good check of 
our definitions and of the numerics of our computation. We have performed this comparison 
and found perfect agreement between the numeric calculation of the real (imaginary) part 
of the correlator and the result obtained via Kramers Kronig relations from the imaginary 
(real) part. 

The DC limit of the above expression gives 



N f N c T 1 (1- ^(1)2)3 | d 2 NfNcT Bd 

axx ~ 1 + & + (l + B2)2' x v 4?r 1 + j g 2 > (5.8) 



which was previously obtained in |30[ and matches the macroscopic calculation in |29[ . We 
have found that there is again a perfect fit. It should be noted that the if) independence 
in the DC Hall conductivity, explicit in the analytic result, is mirrored in the numerical 
solution. This is another very good check that both the conventions for defining the Green's 
function and the numerics themselves are correct. 

In figure |3] we plot the real part of the x — x Green's function for specific values of 
the magnetic field, baryon number density and quark mass as a function of the frequency. 
There is a clear feature appearing in these plots which is not present in the case of zero 
baryon density and/or zero magnetic field. This is the appearance of double peaks, which 
are clearly almost coincident, and the distance apart of which is related to the size of 
B and d. We can study the exact position of these peaks by studying the poles in the 
spectral function in the complex frequency plane, and see clearly the appearance of two 
peaks, emerging from a single peak for B = or d = 0. This effect of the splitting of a 
single quasinormal mode into two is the finite temperature analogue of the Zeeman effect. 
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The breaking of parity symmetry in our system by the magnetic field causes a splitting in 
parity pairs of the quasinormal spectrum. This splitting can be seen explicitly in [[H], |32|] . 

Having studied the spectral function itself we can now look at 'global' properties of 
the Green's functions, defined through sum rules. The first quantity to discuss is the Hall 
angle, given by the ratio of the Hall conductivity to the optical conductivity 

t„(u;)=tanM W ) = ^44- (5 ' 9) 

An example of this quantity is plotted in figure |4j. It has a number properties which differ 
greatly from one material to another. In particular the hall angle tells us about the ratio 



of Hall to drift currents under a given electric field (2.21). 



From the hall angle we can define the hall frequency, ujh, given by ( p. 19 ). In figure || 



we present the numeric results we obtained for the value of the Hall frequency as a function 
of the external magnetic field, the baryon density and the mass of the fundamental degrees 
of freedom. In the examples we present we can see that, for a given B and d, the Hall 
frequency has a maximum value for an specific mass m q . After this maximum is reached, 
the Hall frequency decays with an inverse power of the quark mass. We couldn't find any 
analytic expression for this power or the position of the maximum. 



The strange feature observed in figure |5|(d) is reminiscent of a phase transition [31 



35, 36 1 . We can see that the Hall frequency is sensitive to this change of phase, which 
occurs for m q ~ 0.27 for B = 5 and d = 0.1. It is not a surprise that the Hall frequency 
has a jump at the value of the critical mass, since the system goes to a chiral setup when 
the magnetic field is increased, and we expect the diffusive processes of the plasma to be 
modified. 

6. Conclusions and perspectives 

It would be interesting to find general sum rules for the retarded correlators in the context 
of AdS/CFT. The biggest problem is the evaluation of the asymptotic value G^. As 
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Figure 5: Hall frequency obtained from the Hall angle sum rule as a function of the mass for 
d = 0.1 (blue), d — 0.5 (red) and d = 1 (yellow) and B = 0.5 (top-left figure), B = 2 (top-right 
figure), B — 5 (bottom-left figure) . The Hall frequency tells about the relative strength of the optical 
to Hall conductivities in the large frequency domain and thus we see that there is a, for a fixed B 
and D, a value of the mass for which the relative strengths are closest. Below and above this value 
the optical conductivity dominates the transport of charge carriers in the material. The bottom-right 
figure corresponds to the blue curve in the bottom-left one. We observe the multivaluedness signaling 
the presence of a phase transition for these values of the magnetic field and charge density. 



pointed out in ||, these values are related to expectation values of local operators via 
OPEs in the field theory side. Also, it is expected that the results at weak and strong 
coupling differ considerably, as shown in ||] for J\f = 4 SYM. Comparison with results 
from lattice in the study of thermal plasmas will certainly shed light on whether we should 
trust holographic estimates of transport coefficients from the spectral functions, although 
it is known that the calculation of real time thermal properties on the lattice is extremely 
difficult. In [^] the authors are able to find analytic expressions for a number of sum rules 
in a gravitational context, and it would be interesting to check whether the calculations 
presented here can be merged with their approach. 

This tool may also be useful in the AdS/CMT approach, in which the determination of 
the properties of the AC conductivity is playing an important role. It may also be helpful 
for determining the strength of the delta functions appearing at zero frequency, and in this 
form give information about the gap frequency of the models under consideration. 

In this paper we have calculated sum rules related to both the optical and hall con- 
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ductivities for probe flavors in the plasma of adjoint matter. The results that we obtain 
tell us about how fundamental matter affects these transport properties in the quenched 
approximation. It is know that the probe limit, which has been assumed throughout, does 
not allow one to follow the dynamics of such processes for infinitely long times, as the 
effects of the fundamental matter on the adjoint degrees of freedom will eventually become 
non-trivial. In the preceding timescale however we are able to calculate finite conductivities 
and thus finite plasma frequencies and hall conductivities, both in the DC and AC regimes 
(indeed, as shown, the limit of our AC results agree with the previously calculated DC 
results which are found through classical D-brane solutions in the absence of fluctuations 
and do not rely on the regularization of Green's functions). The fact that the sum rules 
link the well established DC results, to the global properties of the AC results is a strong 
indication that these calculations are correct. 

We have calculated both the plasma frequency and hall frequency which in such a setup 
tell us about the response of fundamental matter to microscopic AC electric currents in 
addition to finite magnetic fields and baryon densities. Although clearly such fine properties 
are not accessible to heavy ion experiments, the techniques outlined here allow for similar 
analysis in lower dimensional systems which should be comparable to condensed matter 
systems. The full analysis of these systems is in progress. 

The fact that we are able to extract highly non-trivial properties of strongly coupled 
systems via holographic sum rules takes us another step closer to calculating measurable 
quantities in real-world systems. The backgrounds discussed in this paper are the simplest 
we could study, but the results here can be extended quite easily to more realistic systems 
which one may eventually be able to engineer in the laboratory. 
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